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Abstract. We present an improved orbit determination 
for the visual pair in the M-dwarf triple system Gliese 866 
that is based on new speckle-interferometric and HST ob- 
servations. The system mass is M = 0.34 ± 0.03 M©. 
The masses of the components derived using the mass- 
luminosity relation are consistent with this mass sum. All 
three components of Gliese 866 seem to have masses not 
far from the hydrogen burning mass limit. 
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1. Introduction 

It is now a well established fact that there exists a large 
number of substellar object s (see e. g. the recent review ar- 
ticle by Oppenheimer et al. 2000). This increases the need 
for a better understanding of stellar properties at the lower 
end of the main sequence. Nearby low-mass stellar systems 
are particularly important for this purpose, because their 
orbital motion allows dynamical mass determinations. 
Here we will consider the nearby triple system Gliese 866 
(Other designations: LHS 68, WDS 2238 5-1519). Using 
speckle interferometry, Leinert et al. ( 1986| ) and also Mc- 
Carthy et al. ( 1987 ) discovered a companion (henceforth 
Gliese 866 B) located about 0'.'4 away from the main com- 
ponent Gliese 866 A. Following observations proved the 
possibility to cover the whole orbit of this binary system 
within a few years by speckle-interferometric observations. 
Based on 16 data points, Leinert et al. ( 1990 , hereafter 
L90) presented a first determination of orbital parame- 
ters and masses. They derived the combined mass of the 
system to be 0.38 ± 0.03 Mq. This was inconsistent with 
values of Ma « 0.14 M© and Mb » 0.11 M© obtained 
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from empirical mass-luminosity relations and stellar inte- 
rior models (see L90 and references therein). 
In the meantime an additional spectroscopic companion 
(hereafter called Gliese 866 a) to Gliese 866 A has been de- 
tected (Delfosse et al. 1999). This is a plausible reason for 
the mentioned mass excess in Gliese 866, but - given the 
fact that only 0.38 Mq now had to be distributed among 
three stars - raises the question for a substellar compo- 
nent in the Gliese 866 system. 

To improve on the mass determination of the compo- 
nents of Gliese 866 we have taken 20 more speckle- 
interferometric observations and one additional HST de- 
termination of relative position. These show the orbit 
of the wide pair: B with respect to Aa. We present an 
overview of the observations and the data reduction pro- 
cess in Sect.g. The results are given in Sect.^, discussed 
in Sect.U and are summarized in Sect.|5|. 

2. Observations and data reduction 

A list of the new observations and their results is given 
in Tablej^. The observations numbered 1, 3, 5 and 7 used 
one-dimensional speckle-interferometry. Details of obser- 
vational techniques and data reduction for this method 
are described in Leinert & Haas (1989). 
All other speckle observations were done using two- 
dimensional infrared array cameras. Sequences of typically 
1000 images with exposure times of « 0.1 sec were taken 
for Gliese 866 and a nearby reference star. After back- 
ground subtraction, flatfielding and badpixel correction 
these data cubes are Fourier-transformed. 
We determine the modulus of the complex visibility (i. e. 
the Fourier transform of the object brightness distribu- 
tion) from power spectrum analysis. The phase is re- 
cursively reconstructed using two different methods: The 
Knox-Thompson algorithm (Knox & Thompson 1974) and 
the bispectrum analysis (Lohmann et al. 1983). Modulus 
and phase are characteristic strip patterns for a binary. 
As an example we show them in Fig. Il| for the observation 
done at 27 September 1996 on Calar Alto. By fitting a bi- 
nary model to the complex visibility we derive the binary 
parameters: position angle, projected separation and flux 
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Table 1. Results of new observations of Gliese 866. Position angle and projected separation are those of Gliese 866 B 
relative to Gliese 866 Aa. The abbreviations in column "Telescope" denote Calar Alto (CA), Kitt Peak National 
Observatory (KPNO), European Southern Observatory (ESO) and New Technology Telescope (NTT). The latter has 
an aperture of 3.5 m. 



No. Date 



Epoch 



Position 
angle [°] 



Projected 
separation 

fmasl 



Filter 



Ib/Iah 



Telescope 
(see caption) 



Camera 



1 


01.09.1990 


1990.6681 


348.0 ± 2.2 


481 ± 18 


K 


0.58 ± 0.03 


CA 3.5 m 


ID 


2 


04.12.1990 


1990.9254 


339.8 ± 2.4 


463 ± 19 


K 




KPNO 3.8 m 


2D 


3 


19.09.1991 


1991.7173 


186.5 ± 3.4 


186 ± 11 


K 


0.55 ± 0.02 


CA 3.5 m 


ID 


4 


29.10.1991 


1991.8268 


155.5 ± 3.8 


200 ± 13 


K 


0.51 ± 0.03 


CA 3.5 m 


1 - 5 jj.m 


5 


16.05.1992 


1992.3751 


15.3 ± 0.5 


244 ± 6 


K 


0.57 ± 0.05 


ESO 3.6 m 


ID 


6 


14.10.1992 


1992.7885 


350.8 ± 1.5 


449 ± 11 


K 


0.58 ± 0.01 


CA 3.5 m 


1 - 5 fim 


7 


10.01.1993 


1993.0273 


348.6 ± 1.4 


480 ± 12 


K 


0.55 ± 0.03 


CA 3.5 m 


ID 


8 


29.07.1993 


1993.5750 


319.3 ± 1.2 


259 ± 11 


K 


0.618 ± 0.071 


ESO NTT 


SHARP 


9 


02.10.1993 


1993.7529 


277.7 ± 0.5 


105 ± 4 


K 


0.60 ± 0.03 


CA 3.5 m 


MAGIC 


10 


01.05.1994 


1994.3313 


74.9 ± 2.7 


126 ± 4 


K 


0.562 ± 0.031 


ESO NTT 


SHARP 


11 


14.09.1994 


1994.7036 


5.9 ± 0.4 


327 ± 4 


K 


0.552 ± 0.009 


CA 3.5 m 


MAGIC 


12 


12.12.1994 


1994.9473 


354.0 ± 0.2 


434 ± 4 


K 


0.551 ± 0.017 


CA 3.5 m 


MAGIC 


13 


13.12.1994 


1994.9501 


354.3 ± 0.3 


439 ±4 


917 nm 


0.603 ± 0.009 


CA3.5m 


MAGIC 


14 


09.07.1995 


1995.5201 


337.3 ± 0.3 


439 ±4 


K 


0.581 ± 0.014 


ESO NTT 


SHARP 


15 


16.08.1995 


1995.6243 


333.1 ± 0.3 


392 ± 4 


K 


0.587 ± 0.022 


ESO 3.6 m 


SHARP 2 


16 


08.06.1996 


1996.4381 


131.6 ± 0.2 


156 ± 1 


583 nm 


0.69 ± 0.06 


HST 


FGS3 


17 


27.09.1996 


1996.7419 


28.0 ± 0.5 


188 ± 4 


K 


0.568 ± 0.022 


CA 3.5 m 


MAGIC 


18 


25.08.1997 


1997.6489 


340.2 ± 0.2 


485 ± 4 


K 


0.593 ± 0.006 


ESO 3.6 m 


SHARP 2 


19 


16.11.1997 


1997.8761 


333.3 ± 0.1 


391 ±4 


K 


0.58 ± 0.01 


CA 3.5 m 


MAGIC 


20 


07.05.1998 


1998.3477 


214.8 ± 0.3 


104 ± 7 


K 


0.558 ± 0.004 


ESO NTT 


SHARP 


21 


10.10.1998 


1998.7748 


101.2 ± 0.2 


133 ±4 


K 


0.497 ± 0.026 


CA 3.5 m 


OMEGA Cass 




Fig. 1. Modulus (left) and 

bispectrum-phase (right) of the 
complex visibility for Gliese 866. 
They have been computed from 
speckle-interferometric data ob- 
tained on 27 September 1996 at the 
3.5m-telescope on Calar Alto in the 
K-band using the NIK array camera 
MAGIC. The bispectrum-phase 
has only been calculated inside 
a circle that corresponds to the 
Nyquist spatial frequency which is 
7.0arcsec~^ for the adopted pixel 
scale. 



ratio. 

To obtain a highly precise relative astrometry which is 
crucial for orbit determination one has to provide a good 
calibration of pixel scale and detector orientation. For the 
speckle observations since 9 July 1995 this calibration has 
been done using astrometric fits to images of the Trapez- 
ium cluster, where precise astrometry has been given by 
McCaughrean & Stauffer (1994). During the previous ob- 
serving runs binary stars with well known orbits were ob- 
served for calibrating pixel scale and detector orientation. 
By doing subsequent observations of these systems and 



calibrating them with the Trapezium cluster we have put 
all speckle observations since July 1993 in a consistent 
system of pixel scale and detector orientation. 



3. Results 

After combining the relative astrometry from L90 and 
our new observations, there are now 37 independent 
data points for the orbital motion of the visual pair in 
Gliese 866. They are plotted in Fig.||, together with the 
result of an orbital fit that used the method of Thiele 
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Table 3. Derived system mass for different subsets of the 
data 
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Y : R. A. [arcsec] 

Fig. 2. biiproved relative orbit of Gliese 866 AB (solid) 
compared to the orbit derived by L90 (dotted). North is 
up and east is to the left. The individual measurements 
are connected to the predicted positions by lines. For the 
two measurements represented by open circles only one co- 
ordinate could be derived from the observation, the other 
one is taken from the predicted orbit. The triangle denotes 
the HST measurement (No. 16 in Table|l|). The periastron 
is indicated by a vertical bar. The line of nodes also is 
plotted. 

Table 2. Orbital elements for Gliese 866 



Orbital element 



Result 



Node 

Longitude of periastron 

Inclination 

Semi majoraxis 

Period 

Eccentricity 

Periastron time 



Q = 162.1° ±0.4° 

uj = -17.7° ±1.1° 

i = 112.4° ±0.5° 

A = (0.346 ± 0.004) arcsec 

P = (2.2506 ± 0.0033) yr 

e = 0.437 ± 0.007 

tp = 1987.236 ± 0.014 



and van den Bos, including iterative differential correc- 
tions (Heintz 197S). The orbital elements resulting from a 



fit to the full data set are given in Table g. Since we don't 
know which node is ascending and which one is descend- 
ing, we choose - as usually is done - fi to be between 
0° and 180°. uj is the angle between the adopted fi and 
the periastron (positive in the direction of motion), and 
i > 90° means clockwise motion. 
Gliese 866 has a trigonometric parallax tt — 289.5±4.4 mas 



van Altena et al. 1995). In the following calculations we 



use the external error of the parallax: 6.8 mas instead of 
the (internal) value given by van Altena et al. This yields 
a distance of 3.45 ± 0.08 pc, a semi majoraxis of 1.19 AU 



Data set 


N 


Msys/MQ 


all data 


37 


0.336 ± 0.026 


only our speckle data 


28 


0.336 ± 0.032 


only our 2D speckle data with 


13 


0.340 ± 0.035 


direct astrometric calibration 







and finally - using Kepler's third law - a system mass 
A/sys = 0.336 ± 0.026 M0. This result remains unchanged 
within the uncertainties if the calculation is done with 
natural subsets of the data (see Table ||). There is partic- 
ularly no significant difference if only the 2D data points 
with good astrometric calibration (see Sect.||) are used. 
Because it covers the longest time span, we take the re- 
sult for the full data set as best values for orbit and system 



4. Discussion 

Most of the uncertainty in system mass is from the par- 
allax error. Further improvements in determining this pa- 
rameter will improve the accuracy of system mass consid- 
erably beyond the present ±7.5 %. 

We also want to get a first estimate of the components' 
masses in order to judge the probability for a substellar 
object in the system. This cannot be done empirically from 
our data, because there are no published radial velocities 
for Gliese 866. Instead we use t he m ass-luminosity relation 
given by Henry & McCarthy (1993) as 



log ^±0.067 = 



-0.1668Mk± 0.5395 



(1) 



for absolute K magnitudes 9.81 < Mk < 7.70. 
This approach is only valid if there is no large shift be- 
tween the spectral energy distributions of the compo- 
nents. To check this we consider the observations taken 
at other wavelengths. At 917 nm the flux ratio is compa- 
rable to that in the K-band and also to the flux ratios in 
other NIR filters given by L90 (Table 3 therein). This sup- 
ports the conclusion of L90 that both (visual) components 
nearly have the same spectral type and effective tempera- 
ture. At 845 nm L90 have measured a higher fiux ratio of 
Ib/Ikb. = 0.83 ±0.1, but this wavelength hes at the edge 
of a strong TiO absorption feature (L90, Fig. 8 therein), 
so an extrapolation of flux ratios into the visible is not 
straightforward. Furthermore Henry et al. ( |1999 ) have ob- 
served Gliese 866 with the F583W filter of the HST. The 
resultant fiux ratio in V is Ib/Iah = 0.69 ± 0.06 and thus 
again close to the values in the near infrared. The fact that 
the fiux ratio Ib/Iah. is nearly constant over a large range 
of wavelengths indicates that all three components of the 
Gliese 866 system have similar effective temperatures and 
thus similar masses. This idea is further supported by the 
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combined spectrum of Gliese 866 (L90, Fig. 8 therein) that 
shows the deep molecular absorptions of an M5.5 star. 
The apparent magnitude of the system is K — (5.56 ± 
0.02) mag (Leggettl992 and references therein). Com- 
bined with the distance given above the absolute sys- 
tem magnitude is (7.87 ± 0.04) mag. The K-band flux ra- 
tios given in TableH result in a mean value of Ib/Iah — 
0.57 ± 0.01. We take the components of the spectroscopic 
pair Gliese 866 Aa to be equally bright. Because we have 
only given qualitative arguments for this assumption, we 
use an error reflecting this uncertainty: /a/^A = 1.0 ±0.5. 
This yields the components' absolute K magnitudes: 



Mk{A) = AfK(a) = (9.11 ± 0.32) mag 

Mk(S) = (8.98 ±0.05) mag (2) 

The resulting masses from Eq.|| then are: 



Ma = Ma = (0.105 ± O.O21)M0 

Mb = (0.110 ±0.018)Mq. (3) 
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The given uncertainties originate from the error of 
the K magnitudes (Eq.H) and the error of the mass- 
luminosity relation itself (Eq. m . The sum of these masses 
is Msys = 0.320 ± 0.035 Mq and is thus within the 
uncertainties consistent with the dynamical system mass 
Msys,dyn = 0.336 ± 0.026 Mg derived above. 



5. Summary 

From a new determination of the visual orbit we have 
given an improved determination of the system mass in 
Gliese 866. With simplifying assumptions we have given 
estimates for the components' masses using the mass- 
luminosity relation by Henry & McCarthy (1993). The 
sum of the components' masses estimated in this way is 
consistent with the dynamically obtained system mass. We 
conclude that there is no substellar object in the triple sys- 
tem Gliese 866 despite the fact that the total system mass 
is only 0.34 Mg. 
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